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A new scenario for the early era of the Universe is proposed. It corresponds to a smooth 
transition between a de Sitter-like phase and a radiation dominated era. We calculate 
the production of gravitons in this model. 
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1. Motivation 

Inflation is the most consistent model to explain the anisotropies in the cosmic 
microwave background radiation (CMBR), and also the origin of the large scale 
structure (LSS) of the Universe. As the Universe evolves, there is, in addition, 
a production of a stochastic background of primordial gravitational waves (GW), 
originated from the vacuum fluctuations. 

In what follows, we set a phenomenological model for unifying the transition from 
the inflationary era to the subsequent radiation dominated epoch of the Universe. In 
particular, we consider a way to extend the framework of the generalised Chaplygin 
gas (GCG)^ to the primordial evolution of the universe. 

This transition is not well understood and it is important to investigate the 
possible signatures in the power spectrum of the stochastic background of GW 
associated with it, giving a new insight on the inflationary model behind the 
accelerated expansion in the early universe. 

2. A GCG inspired scenario for the first stages of the Universe 

This unification is achieved by considering a perfect fluid whose energy density, p, 
has a dependence on the scale factor, a, of the form^ 

B \ 

where A, B are constants determined by the matter content of the universe, and a 
is a free parameter. The conditions, 1 + a < 0, and A, i? > 0, lead to the expected 
interpolation. The conservation of p implies the following equation of state 

1 4 A 

3. The Infiationary Scalar Field 

The inflationary dynamics of this model can be described by a scalar field, 
rolling down its self-interacting potential, V{(j)). This potential can be obtained 
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by identifying p and P of the fluid with those of the scalar field. On the other hand, 
the number of e-folds of expansion since a given mode k exits the horizon (fc — aH) 
during the infiationary era, at = 0c, until the end of infiation can be approximated 

4. Observational Constraints 

The model can be constrained using CMBR/LSS recent measurements. More 
precisely, with the values^ Ps{kc) = (2.45 ± 0.23) x 10"^ and Us = 1.0 ± 0.1 for 
the power spectrum and spectral index of the scalar perturbations, respectively, 
evaluated at the pivot wave number kc = 0.05 Mpc~^. 

Therefore, the energy scale for inflation can be determined through the 
amplitude of the power spectrum for density perturbations, Ps{k) « if;;^^ , 

k—aH 

where = On the other hand, we can analyse the variation of the scalar to 
tensor ratio, r = p = 16 e, where Pt{k) is the tensorial power spectrum, in terms 
of the spectral index of the scalar perturbations^ = 1 — 6 e + 2ri. 

There are two approaches to evaluate r and n^: one using the number of e-folds, 
Nc, and the other one by means of the comoving wavenumber, k — aH. The results 
are presented on the Us — r parameter space in Fig. 1-a, where we see that the 
parameter space is slightly more red than preferred by the observation (blue area) . 
However, Fig. 1-a also shows that r < 1 which is in agreement with the observation.^ 
The most consistent value of a with the observation is —1.024, corresponding to an 
energy scale for inflation of ^ 10^^ GeV. 



5. Gravitational Wave Spectrum 

The GW spectrum can be calculated using the method of the continuous Bogoliubov 
coefficients, a and The dimensionless relative logarithmic energy of the spectrum 
of the GWs, f^GWi at the present time ryo, is related to* /3: 

where 0^ is proportional to the number of gravitons created and lo is the respective 
angular frequency. The same coefficient can be written in terms of two continuous 
functions of time, X and Y ^ which obey the differential equation* 



with X' — —ikY. As initial conditions for the numerical integration of Eq. ([4]); i.e. 
the initial values for ^(77) and Y{r]), we use the exact analytical solution for a dS 
Universe* because our model is quasi-Sitter in the past. For the late-time evolution 
we use the LCDM model defined through the recent measurements of WMAP5. 

In Fig. 1-b we present the GW spectra resulting from the numerical integrations 
using the number of e-folds to constrain the infiationary scale. As can be noticed, the 
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spectrum changes with a, with a simultaneous vertical displacement and a variation 
in the high frequency region. 

The decrease of the plateau height, which can be observed by the difference 
between the solid lines {Nc = 47) and the dashed lines {Nc — 62) for a constant 
a (a given colour), is easily interpreted as a decrease in the inflationary scale. On 
the other hand, we immediately conclude that the increase in the high frequency 
region, when a moves away from —1, is a function of the parameter a only. 




ns f(Hz) 

Fig. 1. Fig. 1-a (left) presents the resulting ns — r parameter space, using the number of e- folds 
of TVc = 47 in blue, and A^c = 62 in purple, or from the comoving wavenumber, fc, in red. The 
black lines are for a constant o-value, parameterised by an increase in the number of e-folds. Fig. 
1-b (right) shows the obtained gravitational wave spectra. 

6. Conclusions 

The GW spectra reveal a consistent picture corresponding to the described smooth 
transition, and represents a significant signature of this modified generalised 
Chaplygin gas. The strong variation of the high frequency range, is directly related 
with the decrease in the maximum of the potential a" /a, when a approaches — 1. 
This issue shows the strong limits to the maximum frequency accepted in our model 
and which will be within the reach of future gravitational- waves detectors like BBO 
and DECIGO,^ for the Khz range of frequencies. In fact, for the most consistent 
values of the a-parameter, the spectrum shows a frequency as low as in the Hz 
region. 
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